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Abstract. High-resolution UCLES/AAT spectra are presented for nine B-type supergiants in the SMC, chosen on the basis 
that they may show varying amounts of nuclear-synthetically processed material mixed to their surface. These spectra have 
been analysed using a new grid of approximately 12 000 non-LTE line blanketed TLUSTY model atmospheres to estimate 
atmospheric parameters and chemical composition. The abundance estimates for O, Mg and Si are in excellent agreement 
with those deduced from other studies, whilst the low estimate for C may reflect the use of the C II doublet at 4267A. The N 
estimates are approximately an order of magnitude greater than those found in unevolved B-type stars or H II regions but are 
consistent with the other estimates in AB-type supergiants. These results have been combined with results from a unified model 
atmosphere analysis of UVES/VLT spectra of B-type supergiants (Trundle et al. 2004 1 to discuss the evolutionary status of these 
objects. For two stars that are in common with those discussed by Trundle et al., we have undertaken a careful comparison in 
order to try to understand the relative importance of the different uncertainties present in such analyses, including observational 
errors and the use of static or unified models. We find that even for these relatively luminous supergiants, TLtJSTY models yield 
atmospheric parameters and chemical compositions similar to those deduced from the unified code FASTWIND. 
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, 1 . Introduction include those of Lennon et al. fl997V Gies & Lambert ("1992), 

hi . Kudritzki et al. (1999 ), McErlean et al. (1999) and Repolust et 

The spectra of early-type supergiants provide an excellent ^ ^ Fitzpatrick & Bohannan Puis 

method for studymg Local Group galaxies due to their high ^r^.^ ^^^^^^ ^^^^ ^^^^^j^ 

intnnsic luminosities (typically 1000 to 10000 Lq). Yet it et al. ( 2004 2005 ) and Venn i 1999J for the Magellanic System. 

IS this extreme luminosity which also gives rise to the com- ^^^^ ^^'^j ^^^^ ^ ^^^^^^ ^ 

plex nature of supergiants, whose stellar wind and extended ^3^^ Monteverde et al. (WTWi and Urbaneja et al. 

atmospheres drive them away from a plane-parallel, Local ^ for M33 and Kaufer et al. (2004), Urbaneja et al. m 



Thermodynamic Equilibrium (LTE) regime. Developments in ^^^^ M|2003b) for other Local Group galaxies. 

theoretical techniques that incorporate non-LTE effects, line j- r.u . / n . . 1 i -,rv/-vo i 

^ ^ Studies of the progenitor O-type stars (see Bouret et al. 120031 

blanketing, sphericity and mass-loss have lead to a resurgance xj-n- » 1 onm o » 1 oaa/i , j d » • » / 

°' ^ ^ ^ , Hilher et al. .2003. , Heap et al. .2004.) and B-type giants (see 

of interest in these objects (for example see Hubeny 1988 ^ , Hnm i t . 1 i ^nm t u ^ j j 1 

, , ^ . „ . , „„„ „ . ... Korn et al. 20021 Lennon et al. 2003t have also provided rele- 

Hubeny et al. .l998; Santolaya-Rey et al 1997 Hillier & Miller , , '—— ' , .. — f,. , . „ 

■ ■ ■' ' ■' ■' vant observations to aid our understanding of the evolution of 

119981. Recent studies of OB-type Supergiants in the Local ^, „ u • 1 r 1 . 

' ' ■'^ ^ the surface chemical composition of early-type stars. 

Group have had a number of aims including understanding the „. iiot 11 j , r.^-T i -iir n 

. . „ ...... T . Since the studies of Jaschek & Jaschek dl 96711. Walborn 



evolution of massive stars, studying the chemical composition j i-,„_^' 1 j- • 1 ' j • 

^ , , , , ,., • , • , (^1972) and Dufton (1972), the dispersion of observed nitro 

of the host galaxy and calibrating the wind momentum versus 



luminosity relation for distance determinations. These studies 



gen abundances in OBA-type stars selected within a partic- 
ular metallicity environment has stimulated continuing inter- 
Send offprint requests to: P.L. Dufton, e-mail: est in the area of massive stars. These changes in the surface 
p . duf ton@qub .ac.uk chemical composition of hot stars are closely coupled to their 
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Table 1. Observational summary of the nine program stars. 
Spectral types are those assigned by Lennon (1997). Apparent 
V-magnitudes are taken from Massey (2002) where available 
and supplemented by Garmany et al. (1987;^) 



Star 


V 


Spec. 


S/N ratio 


Vr 


FWHM 




[mag] 


Type 




km s ^ 


AV78 


11.05 


B1.5Ia+ 


130 


178 


80 


AV215 


12.69 


BNOIa 


90 


159 


115 


AV242 


12.1lt 


Blla 


110 


142 


100 


AV303 


12.78 


B1.5Iab 


100 


186 


65 


AV374 


13.09 


B2Ib 


65 


125 


70 


AV462 


12.54 


B1.5Ia 


115 


126 


70 


AV472 


12.62 


B2Ia 


90 


138 


70 


AV487 


12.58 


BCOIa 


125 


183 


100 


SK191 


11.86^ 


B1.5Ia 


125 


134 


130 



rotational velocities. Current theoretical models require ini- 
tially high rotational velocities in order to mix nucleosyntheti- 
cally processed material to the surface (Heger & Langer 2000 
Maeder & Meynet 2000 2001 1. In turn these models predict 
that B-type supergiants and giants should still be rotating rel- 
atively rapidly with larger rotational velocities in lower metal- 
licity regions where the loss of angular momentum through the 
strong stellar wind is less significant. However as discussed by 
Howarth et al. J1994> . Lennon et al. ( I2003> and Trundle et al. 
J200 4). this appears inconsistent with the observed widths of 
the metal line spectra. Ryans et al. S2002> have attempted to 
distinguish between the different broadening mechanisms (viz. 
microturbulence, rotation and macroturbulence). For a sample 
of Galactic supergiants, they found that the contribution of rota- 
tion to the line-broadening was relatively small, in turn leading 
to large discrepancies with the theoretical predictions. 

In this paper, we present high resolution observations of a 
sample of SMC supergiants chosen from a preliminary explo- 
ration of B-type supergiants by Dufton et al. J2000t . The se- 
lection of targets was made on the basis that they might show 
different degrees of nucleosynthetic processed material at their 
surface. These have been analysed using a new grid of ap- 
proximately 12000 models computed with the non-LTE codes 
TLUSTY and SYNSPEC (Hubenv 1988. Hubeny & Lanz 1995^, 
Hubeny et al. ll998t . These models do not include any contri- 
bution from a stellar wind and therefore, for two targets, we 
have compared our results with those deduced by Trundle et al. 
(2004) from the unified code FASTWIND (Santolaya-Rey et al. 
[1997' with some updates described in Herrero et al. 2002 and 
Repolust et al. 2QQ4j to investigate the magnitude of both ob- 
servational and theoretical uncertainties. In a companion paper, 
we intend to obtain reliable estimates of the actual projected ro- 
tational velocities of our own and other SMC targets using the 
methods discussed by Ryans et al. (I2002> . These should allow 
a better understanding of the relationship between the chemical 
evolution and rotation for massive stars. 



2. Observations and data reduction 

The high-resolution, spectroscopic data presented here were 
obtained during an observing run with the 3.9-m Anglo- 
Australian Telescope (AAT) from 29 September to 1st October 
1998 inclusively. The University College of London Echelle 
Spectrograph (UCLES) was used with the 31 lines mm^^ 
grating and with a TeK IKxlK CCD, providing complete 
spectral coverage between AA3900^900 A at a FWHM res- 
olution of '^O. 1 A. Conditions were excellent throughout the 
three night run, with stellar exposures being bracketed with 
Cu-Ar arc exposures for wavelength calibration. Observations 
were obtained for nine supergiants and these are summarized 
in Table ^ Listed are the stellar V-magnitude (taken from 
Garmany et al. il987 and Massey 2002 1, spectral types (from 
Lennon [l997> . the signal-to-noise (S/N) ratio at approximately 
4500A, the heliocentric radial velocity (vi) and the typical full- 
width-half-maxima (FWHM) of the metal absorption line spec- 
trum. The latter two quantities were estimated as discussed be- 
low, whilst the S/N ratios were deduced by fitting low order 
polynomials to the continuum. Care was taken to include con- 
tinuum regions from different parts of the blaze profile and 
hence these estimates should be considered as conservative, 
particularly given the significant wavelength overlap between 
adjacent orders. 

The choice of targets was based on the spectral types de- 
duced by Lennon i.l991i and the preliminary non-LTE anal- 
ysis of Dufton et al. J2000> of intermediate dispersion spec- 
troscopy. In particular, we attemped to sample a range of lumi- 
nosity types (ranging from Ia+ to lb) and with varying degrees 
of mixing of nucleosynthetic material to the surface. For the 
latter, we included both BN and BC spectral types and stars 
with different ratios of N II to C II equivalent widths (see Fig. 9 
of Dufton et al.). 

The two dimensional CCD datasets were reduced using 
standard procedures within IRAF' (Tody il986i< . Preliminary 
processing of the CCD frames such as over-scan correc- 
tion, trimming of the data section and flat-fielding were per- 
formed using the CCDRED package (Massey .1997 1, whilst 
cosmic-ray removal, extraction of the stellar spectra, sky- 
subtraction and wavelength calibration were carried out us- 
ing the SPECRED (Massey et al. .1992j and DOECSLIT 
(Willmarth & Barnes 1994 1 packages. 

Equivalent widths (EWs) of lines were measured using the 
STARLINK spectrum analysis program DIPSO (Howarth et 
al. I1994> . Low order polynomials were fitted to the adjacent 
continuum regions for normalisation and Gaussian profiles to 
the metal and non-diffuse helium lines using non-linear least 
square routines. As discussed by Ryans et al. J2003> . the pro- 
files of metal absorption lines in the spectra of early-type super- 
giants are well represented by a Gaussian profiles. Tests using 
a different profile shape showed that this assumption was not 
critical to the fitting procedure. 

In addition the central wavelength and widths of the metal 
lines could be used to determine the stellar radial velocity and 

' IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in Tucson 
|.http://iraf. noao.edu I 
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FWHM, which are summarized in Table Typically 10 rela- 
tively strong, isolated lines were used for the former, whilst the 
FWHM were estimated from the Si III multiplet near 4560A, 
which was well observed in all our spectra. Standard deviations 
were typically ±2 km s^^and ±5 km s^ ^respectively, with the 
values for the FWHM being rounded to the nearest 5 km s^^. 
For the hydrogen lines, the identification of the continuum was 
complicated by the presence of a significant echelle ripple. In 
these cases the continua of the adjacent orders were moved to a 
common wavelength scale, fitted using low order polynomials 
and merged. This merged spectrum, which should represent the 
echelle ripple for the order containing the hydrogen line was 
then used to rectify the relevant order prior to normalisation. 
Profiles were then extracted with the continuum levels being 
defined at ±16 A from the line centre. 

3. Spectral analysis 

3. 1. Model atmosphere calculations 

The analysis is based on grids of non-LTE model atmospheres 
calculated using the codes TLUSTY and SYNSPEC (Hubeny 
[19881 Hubeny & Lanz [T9951 Hubeny et al. [1998). Normally 
when using such codes to analyse an observational dataset, the 
approach is to undertake specific calculations, often in the form 
of grids covering a restricted range of atmospheric parameters. 
Such an approach has the advantage that the calculations are 
tailored to the specific problem. Additionally, non-LTE codes 
often exhibit stability problems, and these are particularly se- 
vere when trace ionic species (e.g. Si II at higher effective tem- 
peratures or Si IV at lower effective temperature) are included. 
By limiting the range of atmospheric parameters considered, 
and by selecting the model ions to be included, these stability 
problems can be ameliorated. However this approach is inef- 
ficient when attempting to analyse either large observational 
datasets or different datasets for similar types of objects. 

We are therefore developing an approach based on the cal- 
culation of large grids covering the range of atmospheric pa- 
rameters appropriate to B-type stars. One fundamental assump- 
tion is that the metal line blanketing is dominated by Fe and we 
return to its validity later in this section. Then the structure of 
the atmosphere is defined by the Fe abundance (or 'metallic- 
ity'), the effective temperature (Tcff), gravity (g) and microtur- 
bulence (^t, which affects the amount of blanketing). Currently 
we are calculating 4 grids with metallicities appropriate to our 
Galaxy (i.e. [-^] = 7.5 dex), the LMC (metallicity reduced by 
0.3 dex), SMC (metallicity reduced by 0.6 dex), and low metal- 
licity regimes (metallicity reduced by 1.1 dex). Due to limi- 
tations on the available CPU power we adoped the 'classical' 
TLUSTY ODE Fe model ions for these calculations. The models 
ions for the light elements that have been explicitly included in 
the non-LTE calculations, together with their origins and from 
where they may be obtained, are summarized in Table [2| These 
models come from two main sources, which are discussed in, 
eg. Allende Prieto et al. ( 2003 ) and Lanz and Hubeny ( 2003 1. 

For each metallicity (i.e. Fe abundance), we consider ap- 
proximately 120 different Tcff-log g points (see Fig.[0, which 
are appropriate to B-type stars ranging from the main-sequence 



Table 2. Source of model ions used in generating the non- 
LTE model atmosphere grids. These model ions are discussed 
in Lanz and Hubeny J20031 designated TLUSTY) and Allende 
Prieto et al. (2003 designated Allende). They are available 
from http://tlusty.gsfc.nasa.gov/ and http://hebe.as.utexas.edu/, 
respectively. 



Ion 


Levels 


Source 


Hi 


16 


TLUSTY 


He I 


24 


TLUSTY 


He II 


14 


TLUSTY 


C II 


39 


TLUSTY 


C III 


23 


TLUSTY 


Nil 


51 


Allende 


NlII 


32 


TLUSTY 


On 


74 


Allende 


Oiii 


29 


TLUSTY 


Mg II 


31 


Allende 


Sin 


46 


Allende 


Sim 


74 


Allende 


Si IV 


23 


TLUSTY 


Sii 


14 


TLUSTY 


Siii 


20 


TLUSTY 


Fell 


35 


TLUSTY 


Fe III 


50 


TLUSTY 


Fe IV 


43 


TLUSTY 



to the Eddington limit. At each Tcff-log g point, we then calcu- 
late 25 models covering 5 values of and 5 sets of light metal 
abundances, which are appropriate to the underlying metallic- 
ity regime. This corresponds to approximately 3 000 models 
per metallicity or 12 000 models in total. 

The analysis of a spectrum is then relatively straightforward 
and fast. One selects an appropriate metallicity grid, and then 
analyses spectral features to determine first the atmospheric pa- 
rameters and then the stellar chemical composition. The pro- 
cess is currently semi-automated, and will be further automated 
in the light of experience gained regarding those areas which 
will continue to require user intervention. Further information 
(on for example the spectral lines incorporated in synthesis cal- 
culations and the light element abundances adopted in different 
grids) is available at http://star. pst.qub.ac.uk/| 

The major difficulty with the approach outlined above is 
to maintain stability over the wide range of atmospheric pa- 
rameters (illustrated in Fig. which in turn leads to differ- 
ent ionic species and indeed atomic processes being important 
in different parts of the grid. Closely associated with this has 
been the problem of quality control for the very large num- 
ber of proposed models. However we have now completed the 
SMC metallicity grid and have tested this against existing cal- 
culations using TLUSTY and DETAIL/SURFACE (see Becker & 
Butler TWTWTW Butler "TW Giddings TW Korn 
et al. 2002, McErlean et al. i?99 for details of atomic data 
and codes). Both the overall agreement and expected differ- 
ences (due e.g. to improved model ions) are highly encour- 
aging. Difficulties remain for a small number of spectral fea- 
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4400 4500 
Wavelength (A) 

SMC TLUSTY grid 



4.50 — + 

4.25 - + 

4.00 — + 

3.75 — + 

3.50 - + 

' 3.25 — + 

3.00 - + 

2.75 — + 

2.50 - + 

2.25 — + 

2.00 — + 

1.75 — + 

12 



n„ (/'OOO K) 

Fig. 1. The bottom panel indicates the range of effective 
temperatures (Tes) and gravities (log g) covered by our pro- 
posed grids, with the lowest gravity models being close to the 
Eddington limit. At each point, 100 models are calculated cov- 
ering a range of microturbulence, metallicity and light metal 
abundances. In total, this leads to approximately 12000 mod- 
els. The top panel shows the predicted spectrum for a small 
wavelength region from a typical model with an effective tem- 
perature of 20000K, logarthmic gravity of 3.0 dex, microtur- 
bulence of 10 km s^^and default light element abundances for 
the SMC. 



tures (see http://star.pst.qub.ac.uk/for details), but these are rel- 
atively minor and should not affect the usefulness of the grids. 

These models have been used to calculate spectra (see Fig. 
QJ, which in turn provide theoretical hydrogen and helium line 
profiles and equivalent widths for light metals for a range of 
abundances. Note that as we keep the iron abundance fixed 
within any given grid, we have less extensive data for this ele- 
ment. Tests showed that the grid spacings for our atmospheric 
parameters was sufficiently fine to allow reliable interpolation 
of the theoretical hydrogen and helium profiles at intermedi- 
ate values. For the theoretical metal line equivalent widths, 
these can then be accessed via a GUI interface written in IDL, 
which allows the user to interpolate in order to calculate equiv- 
alent widths and/or abundance estimates for approximately 200 
metal fines for any given set of atmospheric parameters. Ryans 
et al. J2003> reported that the increment of 0.4 dex used in 
our grids was fine enough to ensure that no significant errors 
were introduced by the interpolation procedures. Full theoreti- 
cal spectra are also available for any given model. In summary 
these grids allow a user- friendly non-LTE analysis of hydrogen 



and helium line profiles and of the profiles and equivalent width 
of the lines of light elements over a range of iron line blanket- 
ing and atmospheric parameters appropriate to B-type stars in 
our Galaxy and in Local Group galaxies such as the Magellanic 
Clouds. 

This approach is based on several assumptions. Firstly there 
are those impHcit in the use of the TLUSTY/SYNSPEC pack- 
age. Particularly relevant to the current analysis are the as- 
sumptions of a static atmosphere (which precludes the inclu- 
sion of a wind) and of a plane parallel geometry. We will return 
to these assumptions when comparing the results with those 
generated using the unified code fast wind (Santolaya-Rey et 
al. [T997l Herrero et al. l2002t . Secondly although we have in- 
cluded line blanketing due to Fe, we have excluded that due to 
other iron peak elements, due to limits on the available compu- 
tational power. However experience has shown that the bulk of 
the opacity in the regime of interest is due to Fe (Hubeny et al. 
|l998 ) so this simplification is unlikely to be a source of signifi- 
cant error Thirdly we have also implicitly assumed that the at- 
mospheric structure is not affected by the adopted light element 
(such as CNO, Mg, Si) abundances, which are constrained to 
vary in step. To test this assumption we have considered repre- 
sentative atmospheric parameters spanning our ranges of effec- 
tive temperature, gravity, microturbulence and iron abundance. 
For each set of atmospheric parameters, we have then fixed, for 
example, the O abundance but recalculated models allowing the 
other light element abundances to vary. If the assumption that 
the atmospheric structure does not depend on the light element 
abundances is correct, the oxygen spectrum should not vary 
significantly within these calculations. Our tests show that, ex- 
cluding very weak lines, variations in equivalent width are typ- 
ically less than 1% and are always less than 5%. Additionally 
we have checked that our temperature structures for given sets 
of atmospheric parameters (and iron abundance) do not depend 
significantly on the light element abundances. This is illustrated 
in Fig. 121 where it can be seen that the temperature structures 
agree to typically ±50K with the maximum differences being 
less than ± 150K. Hence we believe that the assumption that the 
atmospheric structure depends to first order on the atmospheric 
parameters and iron abundance is acceptable. 



3.2. Estimation of atmospheric parameters 

The atmospheric parameters of our targets were estimated us- 
ing standard techniques, viz. the silicon ionization equilibrium 
for effective temperaure, the profiles of the Balmer lines for 
gravity and the relative strength of absorption lines of an ionic 
species for microturbulence. These techniques have been de- 
scribed by, for example, Kilian J1992> . McErlean et al. J1999> . 
Korn et al. ('200?! and Trundle et al. r2004i and will only be 
briefly discussed here. However it should be noted that the es- 
timation of the atmospheric parameters is an iterative process 
as they are inter-related. Additionally all the results presented 
below are based on the grid with an iron abundance appropriate 
to the SMC (i.e. 0.6 dex less than solar - Anders & Grevesse 
'1989). Tests were undertaken using other grids and these are 
discussed below. 
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Table 3. Atmospheric parameters and abundances of program stars. Also listed are results for the main sequence SMC target, 
AV 304. The number of lines considered is shown in brackets. 



AV78 AV215 AV242 AV303 AV374 AV462 AV472 AV487 SK191 AV304 



Tcff 21250 26500 21500 18000 18500 19000 19000 27000 20500 27000 
\ogg 2.45 3.00 2.50 2.30 2.50 2.40 2.40 2.95 2.45 3.9 

13 13 13 14 14 16 15 11 16 3 



c 


6.91 (1) 




7.06(1) 


7.15 (3) 


7.07 (1) 


7.00 (3) 


7.06 (3) 


7.24(1) 


7.00 (1) 


7.36 (3) 


a 








0.09 




0.07 


0.12 








N 


7.80 (9) 


7.60 (3) 


7.09 (3) 


7.28 (6) 


7.18(5) 


7.42 (8) 


7.51 (9) 


<7.3 


7.53 (4) 


6.55 (1) 


a 


0.06 


0.06 


0.07 


0.05 


0.05 


0.05 


0.03 




0.02 







7.74 (12) 


8.00(11) 


8.11 (19) 


8.34 (22) 


8.26(12) 


8.12(17) 


8.05 (18) 


8.09 (10) 


8.11 (13) 


8.13(42) 


a 


0.08 


0.08 


0.04 


0.03 


0.03 


0.05 


0.03 


0.04 


0.06 




Mg 


6.79(1) 


6.67 (1) 


6.60(1) 


6.69 (1) 


6.62(1) 


6.73(1) 


6.67(1) 


6.79(1) 


6.71 (1) 


6.77 (1) 


Si 


6.90 (4) 


6.87 (4) 


6.65 (4) 


7.01 (4) 


7.01 (4) 


6.85 (6) 


6.81 (3) 


6.93 (3) 


6.64 (3) 


6.75 (6) 


a 


0.01 


0.02 


0.01 


0.03 


0.02 


0.01 


0.01 


0.03 


0.01 





3.2.1. Effective temperature, T^s 

Initial estimates of the surface gravity were based on the cal- 
ibration of McErlean et al. t,1999j . These were then used to 
estimate the effective temperature normally using the Si III to 
Si IV ionization equilbrium. For AV 374, an additional estimate 
was available from the Si II to Si III ionization equilbrium, with 
the two values agreeing to within 200 K. The Si III multiplet at 
approximately 4560A contains lines with a range of equivalent 
widths and hence the microturbulence could be simultaneously 
constrained (see Sect. I3.2.3> . The final effective temperature 
estimates are summarized in Table |3] rounded to the nearest 
500 K. The quality of the observational data implies that these 
will have an uncertainty of typically ±1 000-2 000 K (see Sect. 
13, while the assumptions implicit to TLUSTY will contribute 
additional uncertainties. We note that Lee et al. ((2004) have 
found that using grids with different Fe abundances changed 
the estimates by typically 500K or less and hence should not 
be a major source of error. 

3.2.2. Surface gravity, g 

The surface gravity was deduced by fitting theoretical profiles 
to the Balmer series features, H/3 to He. As expected the lower 
series members (and in particular H/3) showed evidence of their 
profiles being affected by a stellar wind and in these cases more 
weight was given to the higher series members. The quaUty of 
the fit is shown in Fig.|3]for the target AV472, where the the- 
oretical profiles have been convolved with a Gaussian profile 
to allow for rotational and macroturbulent broadening. Use of 
grids with different Fe abundances did not affect the gravity 
determinations, with the major source of error (of the order of 
±0.2 dex) arising from observational and fitting uncertainties. 
The final adopted gravities, g (in units of cm s^^) are listed in 
TableE). 



3.2.3. IVIicroturbulence, 

The microturbulent velocity was initially determined from 
the Si III multiplet at approximately 4560A, by minimis- 
ing the spread of abundances deduced from its 3 members. 
Subsequently the same procedure was undertaken using the 
rich Oil spectra in our targets. As has been found previ- 
ously (see, for example, McErlean et al. 119991 Vrancken et 
al. liOOOI and Trundle et al. 2004 1, the Oil spectra gave esti- 
mates which were systematically larger by between 5 and 10 
km s^^. Typically 10-20 Oil lines are included but the anal- 
ysis is complicated by them arising from different multiplets. 
Uncertainties in the atomic data and in the magnitude of the 
non-UTE effects may affect such an analysis and hence the 
source of these differences is unclear. We will follow previ- 
ous authors in adopting the results from the silicon multiplet 
but note that the estimates must be considered uncertain by at 
least 5 km s~^(obviously this uncertainty will have a larger ef- 
fect on the abundances derived from strong lines as the weak 
lines are unaffected by the microturbulence, as is apparent from 
Table |3- Tests using grids with different base metallicities 
yielded, as would be expected, effectively identical estimates 
for the microturbulence and the adopted values are again listed 
in TableOl 

3.3. Photospheric Abundances 

The adopted atmospheric parameters (listed in Table |3} were 
used to derive absolute non-LTE abudances for the programme 
stars. Tests were undertaken by comparing the observed equiv- 
alent widths with theoretical values from different metallicity 
grids. The adopted iron abundance was not found to signifi- 
cantly affect the light element abundance estimates (as had also 
been found by Lee et al. "2004) and the results presented in 
Tableware for the SMC grid (metallicity -0.6 dex compared 
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Fig. 3. Observed and theoretical Balmer line proflies for AV472. The theoretical profiles have been convolved with a Gaussian 
profile using the FWHM listed in Table[l] Note the filling in of the Balmer line cores due to wind emission and that this is more 
significant in the lower series members. 



with the solar estimate of Anders & Grevesse ll989> . For ex- 
ample, we re-analysed the metal fines equivalent widths for 
AV215 (one of the targets considered further in Sect.|^ using 
the grids of models with metallicities appropriate to the LMC 
(-0.3 dex) and low metallicity regimes (-1.1 dex). The average 
change in the light element abundance estimates was 0.04 dex, 
with all changes being less than 0. 1 dex. 

Also fisted in this Table are error estimates due to the scatter 
in the abundances deduced from individual features (note that 
there will be other sources of uncertainty including errors in 
the adopted atmospheric parameters or in the physical models 
which we address below and in Sect.|4]respectively). These es- 
timates assume that the errors are normally distributed and are 
simply the sample standard deviations divided by the square 
root of the number of spectral features included. For ions, such 
as Mg II, where only one spectral feature was observed, the 
sample standard deviations from other ions imply that an er- 
ror estimate of 0. 1 to 0.2 dex is probably appropriate. 

VLT observations of the SMC main sequence B-type star, 
AV304 (Rolleston et al. ll996> have been re-analysed with our 
TLUSTY grids (Hunter et al."2004'). This star is particularly suit- 
able for analysis as its spectrum has relatively sharp absorp- 
tion lines. In Table [S] we fist the atmospheric parameter and 
abundance estimates for AV 304 and note that these are in rea- 
sonable agreement with those deduced from the LTE analysis 
of Rolleston et al. J1996> . particularly when the non-LTE cor- 
rections of Lennon et al. (2003 ) are included. We believe that 
these results provide a baseline for the current chemical com- 



position of the SMC and given the methods used are particu- 
larly appropriate for comparison with our results for the SMC 
supergiants. Hence we have undertaken a differential analysis 
for each supergiant relative to AV 304 and these results are pre- 
sented in Table|3 together with error estimates calculated using 
the same methodology as for the absolute abundances. We ex- 
pect that these differential abundances are more reliable as they 
would be less prone to uncertainties in the atomic data. Some 
indirect evidence for this is provided by the smaller error esti- 
mates that are found for the differential, rather than the absolute 
abundance estimates. 

There will be other potential sources of uncertainty in the 
abundance estimates presented in Tables |3l and |4] For exam- 
ple errors in the adopted atmospheric parameters would sys- 
tematically affect the estimates from any given ionic species. 
Therefore in Table |5] we list the changes in the abundance es- 
timates that would arise from increases in the values of T^g, 
log g and by lOOOK, 0.2 dex and 5km s^^ respectively. Note 
that these should be considered as indicative as they may vary 
from line to line in any given ionic species. Values are listed 
for two stars, AV 303 and AV 487, which were chosen as they 
were the coolest and hottest stars in our sample. It is encourag- 
ing to note that given the uncertainties in estimating the mic- 
toturbulent velocities, the errors associated with this quantity 
are relatively small. By contrast the errors associated with the 
effective temperature and gravity are typically 0.1-0.2 dex but 
can be as large as 0.3 dex. These estimates are appropriate to 
the absolute abundance estimates listed in Table |3] In the case 
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Table 4. Differential abundances of program stars relative to the main sequence SMC target, AV 304. The number of lines 
considered is shown in brackets. 

AV78 AV215 AV242 AV303 AV374 AV462 AV472 AV487 SK191 

C -0.38 (1) - -0.23 (1) -0.13 (3) -0.22 (1) -0.36 (3) -0.30(3) -0.05 (1) -0.29(1) 

- - - 0.07 - 0.07 0.10 

N 1.37 (1) 1.00(1) 0.41 (1) 0.79 (1) 0.46(1) 0.93 (1) 0.96(1) <0.7 (1) 0.98(1) 

O -0.43 (12) -0.13 (11) -0.03 (19) 0.20 (22) 0.10(12) 0.01 (17) -0.11 (18) 0.00(10) 0.03 (13) 

a 0.08 0.04 0.03 0.03 0.04 0.03 0.03 0.03 0.03 

Mg 0.02(1) -0.10(1) -0.17(1) -0.08(1) -0.15 (1) -0.04(1) -0.10(1) 0.02(1) -0.06(1) 

Si 0.12(4) 0.08(4) -0.14(4) 0.22(4) 0.21 (4) 0.08(4) 0.04(3) 0.14(3) -0.15(3) 

a(mean) 0.05 0.05 0.05 0.08 0.07 0.04 0.04 0.05 0.06 
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Fig. 2. Differences in the temperature scales as a function of 
the Rosseland mean optical depth for models with an effec- 
tive temperature of 20 000 K, logarithmic gravity of 2.5 dex, 
microturbulence of 10km s^^ and an Fe abundance appropri- 
ate to the SMC. There are five different models with baseline 
SMC element abundances and with light element abundances 
increased and decreased by 0.4 and 0.8 dex. The four plots 
show the difference in the temperature structure of the latter 
compared with that of the baseline abundance model. The small 
differences imply that the temperature structures are effectively 
independent of the light element abundances adopted. 



of estimates listed in Table |4] systematic errors in the deriva- 
tion of the atmospheric parameters might lead to smaller errors 
in the differential abundances and hence these estimates should 
be considered as conservative. 



4. Comparison of FASTWIND and TLUSTY 
analysis tecliniques 



The analysis presented in Sect. |3l is based on static calcula- 
tions using the plane-parallel non-LTE model atmosphere code, 
TLUSTY. Clearly the neglect of the wind must be problematic 
for stars with relatively large luminosities and low gravities 
(see, for example Aller 1956; Kudritzki & Puis 2000 amongst 
others). Indeed features such as the Balmer H/3 line in the spec- 
tra of our current sample show clear evidence of the presence 
of a wind. Hence it is useful to compare our results with those 
deduced from a unified code. In fact VLT/UVES spectra of two 
of our targets (AV 215 and SK 191) have been previously anal- 
ysed by Trundle et al. (2004), facilitating such a comparison. 
The version of the code FASTWIND used by Trundle et al. in- 
cluded some updates from that introduced by Santolaya-Rey et 
al. J 19971 1 for line blanketing and blocking and is briefly de- 



scribed in Herrero et al. (.2002 1 and Repolust et al. ( I2004t . 

Although it is possible to directly compare the results of 
Trundle et al. (hereafter designated FW/VLT) with those found 
here (designated TL/AAT), this has the disadvantage that both 
the observational data and theoretical methods differ between 
the two analyses. To make the comparison easier, we have 
therefore taken the VLT spectra and analysed them using our 
TLUSTY grid and the same procedures as discussed above. This 
analysis is designated TLA'LT and the results for all three anal- 
yses in both stars are summarized in Table|6land are discussed 
below. 

We stress that the comparison discussed below does not at- 
tempt to determine which of the codes has a better descrip- 
tion of the physical processes relevant to the compUcated atmo- 
spheres of B-type supergiants. Instead, it has the more limited 
objective of trying to quantify and understand the differences 
in the estimates for the atmospheric parameters and chemical 
composition obtained for such objects when using codes which 
incorporate different physical assumptions and processes. 
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Table 5. Effects of changes in the adopted atmospheric parameters by ATcfr= 1,000K, Alog g=0.2 and A^t=5km s ^. Results 
are presented for the coolest (AV 303) and hottest (AV 487) stars in our sample. 



Ion 



ATeft- 



Alog g 



ACt 





AV303 


AV487 


AV303 


AV487 


AV303 


AV487 


cn 


+0.10 


+0.21 


+0.01 


-0.18 


-0.01 


-0.01 


Nn 


-0.09 


+0.22 


+0.09 


-0.21 


-0.08 


-0.01 


on 


-0.30 


+0.23 


+0.16 


-0.24 


-0.09 


-0.03 


Mg 11 


+0.14 


+0.14 


-0.07 


-0.17 


-0.01 


-0.07 


Si III 


-0.30 


+0.27 


+0.18 


-0.23 


-0.15 


-0.06 



Table 6. Comparison of atmospheric parameters and abundance estimates for AV215 and SK191 





TL/AAT 


AV215 
TL/VLT 


FWA'LT 


TL/AAT 


SK191 
TL/VLT 


FW/VLT 


Teff 
logg 

?t 


26500 
3.00 
13 


26500 
2.95 
12 


27000 
2.90 
12 


20500 
2.45 
16 


20500 
2.45 
15 


22500 
2.55 
13 


Cii 




7.00 


6.91 


7.00 


6.88 


6.89 


Nil 


7.60 


7.80 


7.96 


7.53 


7.45 


7.63 


On 


8.00 


7.98 


7.97 


8.11 


8.20 


8.20 


Mgll 


6.67 


<7.20 


<7.20 


6.71 


6.87 


6.98 


Sim 
Si IV 


6.87 
6.89 


6.90 
6.85 


7.10 
7.14 


6.64 
6.61 


6.70 
6.51 


6.75 
6.50 



4. /. Observational uncertainties 

A comparison of the TL/AAT and TL/VLT analyses gives an 
insight into the effects of observational uncertainties on the 
resulting model atmosphere analysis. The VLT/UVES spectra 
discussed by Trundle et al. had S/N ratios in the blue spec- 
tra region of 120 (AV215) and 170 (SK 191), compared with 
90 and 125 respectively for the UCLES/AAT data. However 
this comparison is misleading as the VLT/UVES data had been 
re-binned to 0.2A, whilst the pixel size for the UCLES/AAT 
spectra was 0.09A. When the latter is binned to a pixel size of 
O.2A, the S/N ratios of the two spectra become comparable, as 
is illustrated in Fig.|4] 

The agreement between the two analyses is encouraging 
with very similar atmospheric parameters being estimated. For 
example, the effective temperature estimates agree for both 
stars with differences in the gravity and microturbulence es- 
timates being at most 0.05 dex and 1km s^ ^respectively. For 
the abundance estimates, in only three cases is the discrepancy 
greater than 0.1 dex with the maximum discrepancy being 0.2 
dex. For two cases (C 11 and Mg 11 in SK 191) this reflects a 



difference in the equivalent width of the single line used to 
estimate the abundance. For N II, the discrepancy arises both 
from differences in the equivalent width measurements and in 
the sets of lines considered. Hence we conclude that for ob- 
servational data of this quality, the corresponding uncertainties 
are relatively small with abundance estimates having a typical 
accuracy of ±0. 1 dex. 

4.2. TIneoretical uncertainties 

A more complicated but potentially more interesting compari- 
son is that which adopts the same observational data but utilises 
different theoretical approaches, viz. the TL/VLT and FW/VLT 
analyses. The effective temperature estimates differ by 500K 
(for AV 215) and 2000K (for SK 191) with the FASTWIND esti- 
mates being higher We have investigated these differences by 
comparing the temperature structures in the two sets of mod- 
els. Initially the FASTWIND code adopts a TLUSTY temperature 
structure. This is then altered by the considerations of spheric- 
ity and with the requirement that continuity is obtained at the 
transition point between the photosphere and stellar wind. The 
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Fig. 4. Comparison of AAT/UCLES and VLT/UVES spectra 
for the star AV215. The spectra have both been rebinned to 
a pixel size of 0.2A and the VLT/UVES spectrum has been 
shifted by 0.1 in y-direction to aid comparison. The spectral 
region includes the He II line at 4541 A and the Si III multiplet 
ar 4552, 4567 and 4575A 



stronger the stellar wind included in a model at a particular 
temperature the further into the stellar atmosphere this transi- 
tion occurs and once the wind becomes important FASTWIND 
assumes an isothermal atmosphere at a temperature appropriate 
to the given atmospheric parameters. SK191 has a relatively 
strong wind for its spectral type and the requirement for an 
isothermal atmosphere in the wind at a predefined temperature 
minimum leads to the FASTWIND models having lower tem- 
peratures in the atmospheric region where the silicon lines are 
formed (by approximately 1 000 K). Since these Si lines are 
important temperature indicators in this stellar regime, this nat- 
urally affects the effective temperature estimates derived from 
the two codes. For AV215, the wind becomes important at a 
similar Rosseland optical depth to that of SK191. However 
the FASTWIND temperature structure is in close agreement to 
that of the TLUSTY model consistent with the better agreement 
in the effective temperatures estimates for this star. The lat- 



est version of FASTWIND (Urbanei a 120041 Puis et al. I2005t 
incorporates new techniques for delineating the temperature 
structure. These ensure thermal balance for the electrons (see 
Kubat, I'uls & Pauldrach. I1999t . as an alternative to solving 
the radiative equilibrium equation explicitly. Puis et al. show 
that the new temperature structures calculated by FASTWIND 
are in very good agreement with those obtained by the uni- 
fied model atmosphere codes CMFGEN (Hillier & Miller' 1998") 
& WM-BASIC (Pauldrach et al.;2001 ). We have not yet imple- 
mented this version, but it may lead to a more accurate and bet- 
ter understanding of how the wind influences the photospheric 
temperature structure. 

The gravity estimates differ by 0.05 to 0.10 dex, although 
particularly for SK191 this difference will reflect, at least in 
part, the different estimates adopted for the effective temper- 
ature. The microturbulence estimates are in relatively good 
agreement with the maximum difference being 2 km s^^. 
Given the different physical assumptions, model ions and nu- 




Log(T 



Fig. 5. Temperature sti-uctures of TLUSTY and FASTWIND at- 
mospheric models. The solid lines represent the TLUSTY mod- 
els with effective temperatures of 26.5, 22.5 and 20.5 kK as 
marked on the figure. The first represents the model parameters 
for AV 215 and the latter two that for SK 191. The dashed and 
dashed-dotted line represent the equivalent FASTWIND models 
for effective temperatures of 27.0 and 22.5 kK with the appro- 
priate wind parameters for AV215 and Skl91, respectively. All 
other parameters are consistent between these models. Notice 
the Rosseland optical depth at which the wind becomes im- 
portant in the FASTWIND models; identified by the isothermal 
behaviour of the temperature structure. 

merical techniques used in the two sets of calculations, the 
agreement must be considered encouraging. 

It is also possible to compare the abundance estimates but it 
should be noted that this will be complicated by differences in 
the adopted atmospheric parameters. For AV215, the two sets 
of parameters were effectively the same but for SK191, the 
effective temperature estimates differed by 2000K (although 
the differences in temperature structure in the regions where 
the lines are formed may be smaller). Hence we will consider 
the two stars separately. 

For AV215, the C, O and upper limit for Mg abundance 
estimates are in excellent agreement. By contrast, the N abun- 
dance estimates differ by 0.16 dex. This is surprising as previ- 
ous calculations (Becker and Butler . 1989 ) have show that non- 
LTE effects in the N II are relatively small compared with, for 
example, those in Cll (Sigut"1996"l. The differences for Sim 
and Si IV estimates are even larger and range from 0.2 to 0.3 
dex. 

We have attempted to understand the cause of these dif- 
ferences as follows. An increase in the adopted mass-loss rate 
changes the density structure leading to a decrease in the up- 
per photosphere and in the wind. In addition the flux in the 
Paschen continuum is decreased and this causes the reduction 
in the line strength observed in the Si III lines shown in Fig.|6l 
By contrast the Si IV 4116A lines (see Fig.0 show no signif- 
icant change in their line strength, whilst the He II 4541 A line 
strength is actually increased. As the latter represent the pre- 
dominant ionization stage, this suggests a change in the ionisa- 
tion equihbrium, which could be obtained by changing either 
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the temperature or density structures. Test calculations for an 
effective temperature of 27000K show that the difference in 
the temperature structure (of upto 1 OOOK) induced by increas- 
ing the wind from 10^ ''MQyr^^ to 1.35 x lO^^Moyr^^ cannot 
account for the change in the Si line strengths, suggesting that 
it is variations in the density structure which might be more 
important. Without changing the gravity, the only way to alter 
the density structure is to change the /3 parameter which con- 
trols the velocity field in the wind. Increasing/decreasing the 
value of /3 shifts a high mass-loss rate model to higher/lower 
ionisation stages i.e decreasing f3 increases the density. In the 
case of the FASTWIND model for AV215, lowering the /3 value 
to 1.1 (and appropriately adjusting the mass-loss rate so the 
Ha profile is adequately reproduced) decreases the Si IV and 
increase the Si III line strengths, in such a way that the Si abun- 
dance estimate can be reduced by approximately 0. 1 dex. The 
complex dependence of the Si line strengths on the temperature 
structure, microturbulence, and the wind parameters introduces 
some further uncertainty for the silicon abundance estimates 
found from a unified model atmosphere code, and this may be 
reflected in the larger uncertainties quoted by Trundle et al. than 
for the TLUSTY analysis presented here. 

1.05 




0.90 - 



0.85 I I I I 

4530 4540 4550 4560 4570 4580 

Wavelength 

Fig. 6. Effect of stellar- wind on the Si III multiplet at 4552- 
4575A. The spectra are for a FASTWIND model at atmospheric 
parameters appropriate for AV215 (dashed line) and a model 
with the mass-loss rate reduced to 15% of its original value 
(solid line) to simulate a plane-parallel model. Note the af- 
fect on the silicon lines due to the inclusion of the stellar 
wind. These theoretical spectra have been convolved with a 
projected rotational velocity of 91km s~^, consistent with the 
value found for AV215 by Trundle et al. ( l2004b . 

For SK 191, the C, O and Si abundance estimates are in ex- 
cellent agreement, although given the different adopted effec- 
tive temperatures, this may be fortuitous. The Mg values differ 
by 0. 12 dex but this is consistent with the lower effective tem- 
perature used in the TL/VLT analysis. Finally the N estimates 
differ by 0.18 dex with the discrepancy being in the same sense 
as for AV 215. 

We have investigated whether the nitrogen model ions in 
the two calculations could explain this discrepancy. We note 
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Fig. 7. Effect of stellar wind on the Si IV fine at 4416A. The 
spectra are for the same models as shown in Fig. |6l Note that 
for this line, the wind appears to have little effect on the line 
strength. 

that the N III model ion adopted in the TLUSTY grid, spHt the 
gound state term into the two levels with different J values. 
However we have carried out test calculations using a simpler 
N III model ion with a single ground state term and find that this 
has a negligible effect on the predicted N II spectra. Hence we 
do not believe that the adopted N III model ions are the cause 
of the discrepancy. We have also reformatted the TLUSTY N II 
model ion so that it could be used by the FASTWIND code. Test 
calculations with this new model ion at atmospheric conditions 
appropriate for SK 191 and AV215, showed that for the same 
physical conditions, the equivalent widths of singlet N II lines 
(viz. 3995, 4227, 4447 A) were increased by no more than 4% 
whilst the triplets (viz. the 4630 A multiplet) differed by up to 
40%. The latter discrepancy does not have a large effect on the 
mean nitrogen abundance adopted from both analyses but is 
simply reflected in the larger standard deviations for the FAST- 
WIND analysis. 

In the case of SK 191, the discrepancy in the nitrogen abun- 
dance may also reflect the difference in the atmospheric struc- 
tures between the fast wind and tlusty models (see Fig.|5j 
in the region of formation of the nitrogen lines. Additionally 
the N line strengths deduced from FASTWIND models decrease 
as the mass-loss rate is inceased due to changes in the Paschen 
continuum flux. For example reducing the mass-loss rate for 
SK 191 by 15% in test calculations with FASTWIND had the ef- 
fect of increasing the equivalent widths of the nitrogen lines 
by ~ 25%, while for AV215 the omission of the wind in- 
creases the nitrogen equivalent widths by 38%. Hence we be- 
lieve that the lower nitrogen abundances estimates deduced 
from TLUSTY calculations compared with those from the FAST- 
WIND calculation are due at least in part to the inclusion of the 
effect of the stellar wind on the photospheric lines in the latter. 

Given the different physical assumptions, model ions etc, 
used in the two analyses the agreement in the abundance es- 
timates is surprisingly good. For example, the mean of the 
modulus of the differences in the estimates is only 0.10 dex. 
In turn this imphes that even for these luminous la super- 
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giants, an analysis using static atmospheres, although not ideal, 
may still be valid. This is important as in some cases partic- 
ularly for faint extragalactic objects, the observational mate- 
rial to constrain the wind parameters may not be available. 
However our two supergiants are not amongst the most extreme 
la"*" objects, where the effects observed here (such as filling 
in of the lower Balmer lines) might become more significant 
for both the higher Balmer lines and the metal line spectra. 
Additionally for such extreme objects the wind together with 
geometrical effects could influence the temperature structure. 

Although this comparison has been limited to only two su- 
pergiants, it has yielded some important results. Firstly to re- 
liably model the totality of the optical spectra, it is necessary 
to use a model that explicitly includes the stellar wind. This is, 
for example, important for the lower members of the Balmer 
series. However by using features which are not significantly 
affected by the stellar wind (such as higher members of the 
Balmer series) when determining the stellar parameters and 
chemical compositon, an analysis using static models may be 
appropriate. Indeed the complexity of an analysis with FAST- 
WIND or other unified models (as illustrated by our discus- 
sion of the Si spectrum) can introduce additional uncertainties. 
However probably the most important and encouraging result 
from this comparison is that both approaches appear to lead 
to reliable estimates for both the atmospheric parameters and 
chemical compositions. 

5. Discussion 

5.1. Chemical compositions 

In Table0 we summarize the mean abundances obtained from 
our supergiant sample. The errors include the estimated stan- 
dard deviations of these means, which should allow for ran- 
dom errors in the observational data and adopted atmospheric 
parameters. Additionally systematic errors in the effective tem- 
perature and gravity, will to some extent manifest themselves 
as random errors in the abundance estimates due to their dif- 
ferent effects on the abundance estimates in stars with various 
atmospheric parameters (see Table|5|i. By contrast a systematic 
mis-estimation of the microturbulence (see Sect. l3.2.'3l for de- 
tails) would lead to systematic errors in the abundances. We 
have estimated these systematic errors to be of the order of 0. 1 
dex and have included them in quadrature with the random er- 
rors. Note that these error estimates do not include uncertainties 
due to, for example, limitations in the adopted model ions or in 
the physical assumptions adopted in the calculations. Also in- 
cluded in Table0are results for analyses of other targets in the 
SMC, viz. B-type supergiants. Trundle et al. (2004.TLP D); A- 
type supergiants, Venn i fT999ll, V enn & Pryzbilla (I2003a> : NGC 
330 giants, Lennon et al. f2003'); AV304, Hunter et al. (l2004l ) 
and H II regions, Kurt et al. (J999 1. 

Trundle et al. J2004> presented non-LTE analyses of the 
spectra of eight SMC supergiants. Our analysis of nine super- 
giants (two of which are in common with the sample of Trundle 
et al.) effectively doubles the sample size, although our obser- 
vational dataset and methodology limits our analysis to those 
features that can be modelled by a static photosphere. As can 
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Table 8. Mean absolute abundances for the B-type supergiants 
analysed here and by Trundle et al. (2004 TLPD). These values 
represent our best estimates for the mean abundances of B-type 
supergiants in the SMC. 



Element Abundance 



C 7.30±0.04 

N 7.55±0.08 

O 8. 11 ±0.04 

Mg 6.75±0.03 

Si 6.80±0.04 



be seen from TableQ the current results are in good agreement 
with those of Trundle et al. for the elements O, Mg, Si. There 
would appear to be a discrepancy of 0.24 dex in the mean C 
abundances found in the two analyses. However the value of 
Trundle et al. was primarily based on the C II doublet at 4267 A, 
which has been found to give systematically lower abundances 
estimates than other features (Eber & Butler [T988t and leads 
to a mean abundance of 6.96 dex. Trundle et al. then corrected 
this value using the methodology discussed by Lennon et al. 
(2003 ) to obtain the value quoted in Table^] Our C abundance 
estimates are based on both the C II feature at 4267A and the 
lines at 3919 and 3921 A. However if we coiTect our results for 
the former using the same methodolgy, our mean C abundance 
becomes 7.32±0.15 dex, which is in excellent agreement with 
that obtained by Trundle et al. and indeed with other analyses 
of SMC objects. 

For nitrogen, our mean abundance is 0.25 dex lower than 
that of Trundle et al., although it is in reasonable agreement 
with that of other evolved SMC targets. Part of this difference 
might be due to the possible systematic difference of approxi- 
mately 0.15 dex found in the N abundance estimates deduced 
from the TLUSTY and FASTWIND analyses discussed in Sect. 
0] However as discussed by Trundle et al., B-type supergiants 
show a wide range of nitrogen enhancements (compared to un- 
evolved objects) and hence at least part of this difference may 
reflect real variations in the mean nitrogen abundances for the 
two samples. 

In Table |8] we list the mean abundances found for B-type 
supergiants when the current results are combined with those 
of Trundle et al. For elements, N, O, Mg, Si, the estimates have 
been simply averaged, whilst for C, our estimates have been 
corrected as discussed above. The error estimates assume that 
the errors follow a normal distribution (i.e. they are the standard 
deviation of the individual estimates divided by the square root 
of the number of measurements). As such they will not include 
any systematic errors due to, for example, limitations in the 
physical assumptions. Such errors are very difficult to assess 
but we note that in Sect. l4.2l the adoption of different theoretical 
approaches yielded differences in the abundance estimates of 
typically 0.1 dex. These values in Table |8l represent our best 
estimates for the mean abundances of B-type supergiants in the 
SMC. 
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Table 7. Mean absolute abundances for the program stars. Also listed are results for analyses of other targets in the SMC, viz. 
B-type supergiants. Trundle et al. ©04 TLPD); A-type supergiants, Venn ( 1999 Venn & Pryzbilla (2003ai; NGC 330 giants, 
Lennon et al. (2003i; AV304, Hunter et al. (2004i; Hll regions Kurt et al. ( 1999). Note a correction has been applied to the C 
abundance in the TLPD analysis, due to the problematic 4267 Aline. If such a correction was applied to our estimate excellant 
agreement is found (see discussion in Sect.|5}- 



Objects 


AB-type supergiants 


B-type stars 


Hit 


Analysis 


this work 


TLPD 


A-type 


NGC330 


AV304 


regions 


C 


7.06±0.12 


7.30±0.09 




7.26±0.15 


7.36±0.12 


7.53±0.06 


N 


7.42±0.15 


7.67±0.27 


7.52±0.10 


7.51±0.18 


6.55±0.18 


6.59±0.08 


O 


8.09±0.12 


8.15±0.07 


8.14±0.06 


7.98±0.15 


8.12±0.10 


8.05±0.05 


Mg 


6.70±0.10 


6.78±0.16 


6.83±0.08 


6.59±0.14 


6.77±0.16 




Si 


6.85±0.12 


6.74±0.11 


6.92±0.15 


6.58±0.32 


6.74±0.19 


6.70±0.20 



The elements, O, Mg and Si, are unlikely to have been af- 
fected by the mixing of nucleosynthetic material to the surface 
and the agreement between these mean supergiant abundances 
and those for AV 304 is excellent with differences of 0.06 dex 
or less. This provides indirect evidence that the methods used 
for the supergiant spectra are reliable. For C, the supergiants's 
estimate is smaller than that for AV 304 by 0.06 dex, with N be- 
ing typically enhanced by 1.0 dex. As discussed by Trundle et 
al. this is consistent with processed material having been mixed 
to the surface and with the predictions of Maeder & Meynet 
J2001> . Unfortunately Maeder & Meynet only tabulate abun- 
dances ratios, whilst their initial abundances are one fifth so- 
lar, corresponding to a N abundance of 7.3 dex, significantly 
greater than that found in AV 304 and SMC H II regions. Hence 
it is not straightforward to directly compare our absolute abun- 
dances with existing stellar evolutionary calculations. However 
it is interesting to note that stars, such as AV 78, have currently 
N abundances that are larger than the sum of their presumed 
initial C and N abundances. Thus if all the additional nitrogen 
has been formed via hydrogen burning some depletion of oxy- 
gen would also be required. Indeed besides a very high N abun- 
dance, AV 78 also exhibits relatively low C and O abundance 
estimates. 

As discussed by, for example. Trundle et al. i2004> . cur- 
rent stellar evolutionary calculations require large initial stel- 
lar rotational velocities in order to obtain significant N en- 
hancements at the surface of B-type supergiants. Additionally, 
it is predicted that these objects will still be rotating relatively 
quickly; for example, the models of Maeder and Maynet ('2001 1 
with initial rotational velocities of 300km s^^and masses of 20 
to 60 solar masses have predicted rotational velocities of 100 
to 200km s~^for effective temeparture appropriate to early B- 
type supergiants. However Trundle et al. observed metal line 
widths, which would imply projected rotational velocities of 
50-90km s^^and even allowing for projection effects, these 
would appear incompatible with the predictions. We find sim- 
ilar inconsistencies with evolutionary models with the FWHM 
listed in Table^again implying projected rotational velocities 
of less than 100 km s^^. 



The discrepancy may be greater than this simple compar- 
ison implies, as it assumes that the observed line widths are 
dominated by rotation. As discussed by Howarth et al. ( I1997> . 
the lack of any B-type supergiants with narrow lines implies 
that another broadening mechanism, as well as rotation, must 
be present. This was confirmed by Ryans et al. (|2002^, who 
analysed very high quality spectra of Galactic B-type super- 
giants to distinquish between the effects of rotation and tur- 
bulence. They found that turbulence was the dominant mech- 
anism with estimates of projected rotational velocities being 
typically 10-20 km s^^. For the current SMC dataset, it was 
found that Gaussian profiles gave a good fit to the observed 
metal line spectra again implying that rotation was not the dom- 
inant mechanism. Additionally, analyses of the spectra of O- 
type main sequence stars normally indicate N enrichment in 
their atmospheres. For example Heap et al. (2004) found that 
for seventeen O-type stars, fourteen of them exhibited N en- 
richment. Hence it is clear that the mixing process operates at 
an early evolutionary stage and either that it occurs at relatively 
small rotational velocities or that the velocity breaking between 
B-type supergiants and their O-type precursors is greater than 
predicted by current evolutionary models. 

Venn (1999) and Venn and Przybilla (2003a) have esti- 
mated N abundances for SMC A-type supergiants. As can be 
seen from Table their mean value is in excellent agreement 
with the value found by combining our two samples (see Table 
|8}. In turn this implies that little further enrichment occurs as B- 
type supergiants evolve into A-type supergiants, which is con- 
sistent with the predictions of evolutionary calculations (see, 
for example, Maeder and Mavnet l2001t 



6. Conclusions 

We have presented a grid of 12000 non-LTE models cover- 
ing the range of atmospheric parameters appropriate to B-type 
stars. These should allow efficient and reliable analysis of such 
objects in environments with different metallicities. In this pa- 
per, the grid has been used to analyse the spectra of 9 SMC B- 
type supergiants to obtain atmospheric parameters and chem- 
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ical compositions. The principal results of this work can be 
summarized as follows: 

1. The abundance estimates for O, Mg and Si are in excel- 
lent agreement with those deduced from unevolved stars 
and H II regions and provide indirect support for the reli- 
ability of the methods. The high N abundances found in 
other evolved objects (see, for example, Korn et al. 120021 
Trundle et al. 2004 VennlWl are found here and for the 
most extreme cases imply that both the CN and ON cycles 
must have been operating. 

2. For two stars we have compared the results obtained us- 
ing our TLUSTY grid with those deduced using the unified 
code FASTWIND. Even for these luminous la supergiants 
the agreement is excellent with discrepancies being of a 
similar magnitude to the observational uncertainties. 

3. Our estimates of the upper limits for the projected rota- 
tional velocities appear inconsistent with those required by 
current evolutionary models showing significant enhance- 
ments of N at the stellar surfaces. This clearly warrants fur- 
ther investigation into the actual projected rotational veloc- 
ity of B-type SMC supergiants (rather than just upper lim- 
its) and in a companion paper we analyse using the methods 
discussed by Ryans et al. ( 2002 ) both the dataset presented 
here and that of Trundle et al. (2004i . 
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